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ABSTRACT 



Context. Young radio-loud active galactic nuclei form an important tool to investigate the evolution of extragalactic radio sources. To 
study the early phases of expanding radio sources, we have constructed CORALZ, a sample of 25 compact (6 < 2") radio sources 
associated with nearby (z < 0.16) galaxies. 

Aims. In this paper we determine the morphologies, linear sizes, and put first constraints on the lobe expansion speeds of the sources 
in the sample. 

Methods. We observed the radio sources from the CORALZ sample with MERLIN at 1.4 GHz or 1.6 GHz, the EVN at 1.6 GHz, and 
global VLBI at 1.6 GHz and/or 5.0 GHz. 

Results. Radio maps, morphological classifications, and linear sizes are presented for all sources in the CORALZ sample. We have 
determined a first upper limit to the expansion velocity of one of the sources, which is remarkably low compared to the brighter 
GPS sources at higher redshifts, indicating a relation between radio luminosity and expansion speed, in agreement with analytical 
models. In addition we present further strong evidence that the spectral turnovers in GPS and CSS sources are caused by synchrotron 
self-absorption (SSA): the CORALZ sources are significantly offset from the well-known correlation between spectral peak frequency 
and angular size, but this correlation is recovered after correcting for the flux-density dependence, as predicted by SSA theory. 

Key words. Galaxies: active - Radio continuum: galaxies 



1. Introduction 

Young radio-loud active galactic nuclei (AGN) are ideal objects 
to study the trigger of AGN activity, and the early evolution 
of classical double radio sources. Young radio sources can be 
recognised by their small angular (physical) sizes, their inverted 
spectra at low frequencies, and their often compact symmetric 
morphologies, with jet or lobe-like structures on either side of 
a central core. Depending on the observed characteristics, they 
are called Compact Steep Spectrum (CSS) sources (Peacock 
& Wall [19821 Fanti et al. [19901 O' Dea [19981 , Gigahertz 
Peaked Spectrum (GPS) sources (Q'Dea [T998l), and/or C ompac t 
Symmetric Objects (CSO, Wilkinson et al. [1994l Conwav l2D02l) . 

Samples of CSS, GPS, and CSO have been used in a statisti- 
cal way to constrain the early evolution of radio sources. Relative 
number statistics indicated that at a young age, radio sources 
must have significantly higher radio luminosities than at a later 
stage (O'Dea & Baum [l99rjl Readhead et al. [19961 Fanti & Fanti 
2003). This decrease in radio luminosity may be preceded by a 
period of luminosity increase (Snellen et al. 2000b; 20031. 

Very interesting results come from multi-epoch Very Long 
Baseline Interferometry (VLBI) observations of individual GPS 
and CSO. They show, for a handful of powerful objects the 
proper motions of the lobes, from which their dynamical ages 
can be estimated (Owsianik & Conway 119981 Polatidis & 
Conway 2003 and references therein). Separation velocities be- 



tween the opposite extremities of the sources of up to 0.4/z _1 c 
have been measured, corresponding to dynamical ages in the 
range of a few hundred to a few thousand years. 



The CORALZ sample: The most nearby young radio sources 
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Snellen et al. (2004]! have selected a sample of young radio 
sources, CORALZ (COmpact RAdio sources at Low Redshift), 
with the aim of obtaining an unbiased view of young radio-loud 
AGN in the nearby universe. The sources in the CORALZ sam- 
ple are significantly closer than the archetypal CSO and GPS 
sources available in the literature. Proper motions should be 
easier to detect in these nearby sources, since similar intrinsic 
expansion velocities would result in significantly higher angu- 
lar motions compared to previous studies. Moreover, in general 
these radio sources are less powerful than the archetypal CSO 
and GPS sources, which will allow luminosity-dependent com- 
parison studies. 

The CORALZ programme was also initiated to investi- 
gate possible selection biases in GPS and CSO samples, since 
CORALZ is thought to be much less biased, due to its rela- 
tively clean selection criteria. The sample has been selected on 
flux density (Si.4GHz > 100 mJy) and angular size (6* < 2"). 
Using the Very Large Array (VLA) Faint Images of the Radio 
Sky at Twenty-centimeters (FIRST) survey (White et al. 119971 1, 
the optical Automated Plate Measuring machine (APM) cata- 
logue of the first Palomar Observatory Sky Survey (POSS-I) 
(McMahon & Irwin 1991 ), and follow-up observations, all radio 



2 



N. de Vries et al.: VLBI observations of the CORALZ sample 



sources identified with bright galaxies were selected (red mag- 
nitude of e < 16.5 mag or blue magnitude of o < 19.5 mag). 
Originally, four of these 28 sources where excluded from the 
CORALZ sample, since the offset between the radio position 
and its optical counterpart was suspiciously large (A pos > 2"), 
suggesting that the radio source had been erroneously iden- 
tified with a random foreground galaxy. In the mean time, 
based on the Sloan Digital Sky Survey (SDSS) Data Release 
6 (Adelman-McCarthy et al. 2008), one of these identifications 
actually turned out to be correct (CORALZ Jl 15000+552821 
and SDSS Jl 15000.08+552821.4, with an offset of 0.5"), rais- 
ing the number of radio sources with reliable optical identifica- 
tions to 25. These 25 compact radio sources will be referred to as 
the CORALZ sample. The three sources that still show anoma- 
lously large offsets between the radio and optical position, are 
referred to as additional sources. Although the CORALZ sam- 
ple is statistically complete down to the flux density limit, the 
limit corresponds to different redshifts, since the host galaxies 
are selected on optical magnitude and have a range in absolute 
magnitude. However, simulations show that the sample is 95% 
statistically complete in the redshift range 0.005 < z < 0.16 
(Snellen et al. 2004), which, including the new SDSS identifica- 
tion, now contains 18 sources. In the rest of this paper we will 
refer to these 18 sources as the CORALZ core sample. A de- 
tailed description of the selection process and additional radio 
observations can be found in Snellen et al. (120041) . Nearly all 
(17/18) of the sources in the CORALZ core sample are classi- 
fied as GPS and CSS sources, confirming that these are likely 
to be young radio sources. This also confirms that the spectral 
shape of a radio source is a good criterion to select very young 
radio (GPS) sources in general. 

In this paper we present VLBI and MERLIN observations of 
the CORALZ sample. In Sect.|2]the observations and data reduc- 
tion are described; the data analysis and modeling is described in 
Sect. [3] The results are presented and discussed in Sect. |H and 
we conclude with Sect. Throughout the paper we adopt the 
cosmological parameters as found by WMAP5 (Komatsu et al. 
120081 H = 70.1 km s 1 Mpc 1 , Q A = 0.721, Q m = 0.279). 

2. Observations and data reduction 

We observed radio sources from the CORALZ sample using 
a range of arrays at several frequencies; MERLIN, the Multi- 
Element Radio Linked Interferometer Network, at 1.4 GHz or 
1.6 GHz, the European VLBI Network (EVN) at 1.6 GHz, and 
global VLBI at 1.6 GHz and/or 5.0 GHz. The combination of 
telescope array and observing frequency were chosen to have a 
wide range in angular resolution, required to match the expected 
wide range in angular size of the sources in the sample. An initial 
guess of the angular sizes of the individual sources was obtained 
from the position and strength of the radio spectral peak. In this 
way the sample was divided into three sub-samples of expected 
'large', 'intermediate', and 'small' size sources. Details of the 
different observing runs are given below. 

MERLIN observations 

The subsample of 'large' radio sources (seven sources, in- 
cluding four from the CORALZ core sample) was ob- 
served with MERLIN using the following six telescopes: 
Defford, Cambridge, Knockin, Darnhall, Jodrell Bank Mk2, and 
Pickmere. The observing times ranged from 9 to 17 hours, alter- 
nating between the science targets (3 min.) and the radio sources 
used as their phase-references (2 min.), resulting in on-source 



integration times ranging from 5 to 10 hours. The largest pro- 
jected baselines were typically 200 km, resulting in a resolu- 
tion of about 150 mas using uniform weighting. These observa- 
tions were carried out between November 9, 2000, and January 
2, 2001, at either 1.4 GHz or 1.6 GHz. 

EVN observations at 1.6 GHz 

The twelve 'intermediate' size CORALZ sources (includ- 
ing ten sources from the CORALZ core sample) were ob- 
served with the EVN at 1.6 GHz using the Effelsberg, Jodrell 
Bank, Medicina, Noto, Torun, Westerbork, and Onsala anten- 
nas between February 14, 2000, and May 28, 2001, except 
for two sources, CORALZ J134158+541524 and CORALZ 
J150805+342323, which were observed on February 19, 2002, 
using the MERLIN antenna at Cambridge instead of Noto. Each 
source was typically observed for 10 x 13 minutes, spread in 
time to obtain optimal UV coverage. The largest projected base- 
lines were typically 1400 km, resulting in a resolution of about 
20 mas using uniform weighting. J0650+6001 and J1740+5211 
were observed as primary calibrators. The data correlation was 
performed at Socorro on the VLBA correlator. The Astronomical 
Image Processing System (AIPS) has been used for editing, a- 
priori calibration, fringe-fitting, self-calibration and imaging of 
the data. For fringe-fitting, the AIPS task FRING was used with 
a point source model, a solution interval of 6 minutes, and a 
standard signal-to-noise ratio cutoff of 5. Several iterations of 
phase (or amplitude & phase) self-calibration and imaging were 
performed using the AIPS tasks CALIB and IMAGR, until the 
image quality converged, with negligible negative structure and 
a low noise level. 

Global VLBI observations (2000) at 5.0 GHz 

The eight 'small' CORALZ sources (including three sources 
from the CORALZ core sample) were observed on March 2, 
2000, with an 18 station array at a frequency of 5.0 GHz, us- 
ing the EVN antennas Effelsberg, Jodrell Bank, Medicina, Noto, 
Torun, Westerbork, Onsala, and Shanghai and the Very Long 
Baseline Array (VLBA) antennas St. Croix, Hancock, North 
Liberty, Fort Davis, Los Alamos, Pie Town, Kitt Peak, Owens 
Valley, Brewster, and Mauna Kea. Each source was typically 
observed for 3x11 minutes, spread in time to obtain opti- 
mal UV coverage. The largest projected baselines were typically 
9000 km, resulting in a resolution of about 2 mas using uniform 
weighting. J0650+6001 and J1740+5211 were observed as pri- 
mary calibrators. The data correlation was performed at Socorro 
on the VLBA correlator. The data reduction was carried out the 
same way as the EVN observations. 

Global VLBI observations (2004) at 1.6 GHz and 5.0 GHz 

The fourteen most compact sources from the CORALZ core 
sample were observed for a second time on May 24, 2004, and 
June 5, 2004, at 5.0 GHz and/or 1.6 GHz, respectively. A similar 
observing strategy was used as the previous global VLBI obser- 
vations, with observing times of 3 x 1 1 minutes. The largest pro- 
jected baselines were typically 9000 km, resulting in a resolution 
of about 2 mas at 5.0 GHz, and 5 mas at 1.6 GHz, using uniform 
weighting. Both experiments were conducted with global VLBI, 
using the EVN antennas Effelsberg, Jodrell Bank, Medicina, 
Noto, Torun, Westerbork, and Onsala and the VLBA antennas 
St. Croix, Hancock, North Liberty, Fort Davis, Los Alamos, Pie 
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Town, Kitt Peak, Owens Valley, and Brewster. J1740+5211 and 
4C39.25 were observed as primary calibrators. The data correla- 
tion was performed at Joint Institute for VLBI in Europe (JIVE). 
The data reduction was carried out the same way as the EVN 
observations. 

All details of the observations are given in the appendix in 
Table IB. fl in Col. 1 IAU name of the source, in Col. 2 the figure 
in which the map is presented, in Cols. 3-7 the observing date, 
frequency, bandwidth, and array, in Col. 8 the peak flux density, 
in Col. 9 the r.m.s. noise level, and in Cols. 10-12 the major and 
minor axis and the position angle of the restoring beam. 

3. Data analysis 

3. 1 . Source structure and angular sizes 

The observed radio source structures were characterised using 
the AIPS task JMFIT, which fits (by least-squares) up to four 
Gaussian components to an image subsection. A first guess of 
the input parameters (position, flux density, and shape of the 
components) was made using the INPFIT procedure. The results 
for all observations of all sources are given in Table IB.2I with 
in Col. 1 the source IAU name, in Col. 2 the figure in which the 
corresponding map is presented, in Col. 3 the observation date, 
in Col. 4 the observing frequency, in Cols. 5-6 the relative posi- 
tion of the component, in Cols. 7-9 the flux density, deconvolved 
size and position angle of the component. The angular sizes of 
the radio sources were determined by measuring the distance 
between the outer edges of the two outermost components in 
the deconvolved image. In the case of unresolved components, 
the positions of the components were used. For unresolved or 
barely resolved sources the deconvolved major axis as returned 
by JMFIT was used as an estimate of the angular size of the radio 
source. 



J073934+495438) is practically unresolved, and a second 
source (CORALZ J131739+41 1545) shows a complex mor- 
phology, making the analysis of expansion velocities for these 
sources impossible at this stage. The third source (CORALZ 
J083 139+460800) on the other hand exhibits a compact dou- 
ble morphology with two strong, unresolved components sep- 
arated by 4.4 mas. The relative positions of the two compo- 
nents were determined using the AIPS tasks JMFIT as well as 
MAXFIT, which fits a quadratic function to a selected region of 
the map to determine the position of an extremum. It is always 
difficult to obtain error estimates for component positions from 
VLBI observations, which are particularly important in case of 
slow expansion speeds and upper limits. When three or more 
epochs of observations are available, positional errors are often 
derived from the scatter around a linear fit to the component mo- 
tion (Owsianik & Conway Tl 9981) . although in this case no inde- 
pendent goodness-of-fit can be determined. Since we only have 
two epochs of observations, we use the following formula from 
Fomalont ( 1999) to derive estimates of the error in the compo- 
nent position: 



where cr rms is the post-fit r.m.s. error of the map, and d and / pea k 
are the size and the peak intensity of the component. We have 
also used Monte Carlo simulations to check our error estimates, 
by repeatedly fitting the peak position of a gaussian with ran- 
dom noise added. The error estimates from the approximating 
formula from Fomalont ( 1999) and from the Monte Carlo simu- 
lations were very similar. We do note, however, that these error 
estimates only take into account the errors in the analysis of the 
radio maps, not the possible biases introduced during the syn- 
thesis of these maps. 



3.2. Morphological classification 

We also made an attempt to morphologically classify the radio 
sources by eye, using the source structures at different frequen- 
cies. We identify five classes: 

1 . Unresolved (U). The source is unresolved or barely resolved. 

2. Compact Double (CD). The source exhibits two components 
with a similar spectral index and/or flux density. 

3. Compact Symmetric Object (CSO). A source with a range 
of components, with the most flat spectral component not 
located at one of the extremities of the source. 

4. Core-Jet (CJ). A source with two or more components, sig- 
nificantly different in flux density and/or spectral index, with 
the most flat spectral component located at one of the ex- 
tremities of the source. 

5. Complex (CX). A source with a complex morphology, not 
falling in one of the above categories. 

Each identification is followed by a '?' if the classification is 
particularly uncertain. The results are shown in TableQ] and dis- 
cussed in Sect. 14.21 

3.3. Expansion velocity 

For three sources we have obtained two epochs of VLBI ob- 
servations (Fig. IB. 21 ), observed with the same telescope ar- 
ray at the same frequency, for which we could potentially es- 
timate expansion velocities. One of these sources (CORALZ 



4. Results and discussion 

4. 1 . Source structure 

Contour maps of the CORALZ sample at 1.4, 1.6, and/or 5.0 
GHz, using the EVN, global VLBI, and MERLIN arrays, are 
shown in the appendix. Fig. IB. II shows contour maps of those 
eight sources from the CORALZ core sample, observed with the 
EVN at 1.6 GHz in 2000/2001 (left column), and with global 
VLBI in 2004 at 1.6 GHz and 5.0 GHz (middle and right col- 
umn, respectively). Fig. IB. 21 shows contour maps of the three 
most compact sources in the CORALZ core sample, obtained 
with global VLBI at 5.0 GHz in 2000 (top row), and 2004 (bot- 
tom row). Fig. IB.3l shows radio maps of the six more extended 
sources in the CORALZ core sample, taken with MERLIN in 
2000 at 1.4 GHz or 1.6 GHz, or with the EVN at 1.6 GHz. Fig. 
IB. 41 shows contour maps of other sources from the CORALZ 
sample, that are not part of the core sample, obtained with global 
VLBI at 5.0 GHz or with MERLIN at 1.4 GHz or 1.6 GHz. 
These objects are nearby GPS/CSS radio sources, but are not 
included in the CORALZ core sample, because they either have 
slightly too low flux densities (S 1.4 ghz < 100 mJy), or are too 
distant (z > 0.16). Fig. lB.5l shows contour maps of the additional 
sources, which have been excluded from the CORALZ sample 
because the optical and radio positions are significantly offset, 
indicating that the radio sources have been erroneously identified 
with a random foreground galaxy. For three sources (CORALZ 
J073328+560541, CORALZ J103719+433515, and CORALZ 
J160246+524358), the S/N at 5.0 GHz on the longer baselines 
was too low for fringe-fitting, and we used the EVN antennas 
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Table 1. Source characteristics of the CORALZ sample. 



IAU Name 


z 


S 1.4GHz 


Ls.OGHz 


Vpeak 


e 


LLS 


Morph. 






(mJy) 


(W Hz" 1 ) 


(MHz) 


(mas) 


(pc) 




The CORALZ core sample 














J073328+560541 


0.104 


394 


24.68 


460 


47 


90 


CSO 


J073934+495438 


0.054 


107 


23.63 


950 


2* 


2* 


u 


J083 139+460800 


0.127 


131 


24.62 


2200 


9 


20 


CD 


J083637+440109 


0.054 


139 


23.66 


< 150 


1600 


1700 


CSO? 


J090615+463618 


0.085 


314 


24.49 


680 


31 


49 


CSO 


J102618+454229 


0.153 


105 


24.55 


180 


17 


45 


CSO 


J103719+433515 


0.023 


129 


22.96 


< 150 


19 


9 


CSO 


Jl 15000+552821 


0.139 


143 


24.57 


< 230 


41 


100 


u 


J120902+411559 


0.095 


147 


24.26 


370 


20 


35 


CSO 


J131739+411545 


0.066 


249 


24.37 


2300 


4 


5 


cx 


J140051+521606 


0.116 


174 


24.36 


< 150 


150* 


320* 


u 


J140942+360416 


0.148 


143 


24.45 


330 


27 


70 


CJ? 


J143521+505122 


0.099 


141 


24.20 


< 150 


150* 


270* 


u 


J 150805+342323 


0.045 


130 


23.35 


< 230 


170 


150 


CD? 


J160246+524358 


0.106 


576 


24.75 


< 150 


180 


350 


CSO 


J161 148+404020 


0.152 


553 


25.03 


< 150 


1300 


3400 


cx 


J 170330+454047 


0.060 


119 


23.54 


< 150 








J171854+544148 


0.147 


329 


24.86 


480 


68 


175 


CSO? 


Other nearby sources in the CORALZ sample 










J093609+331308 


0.076 


55 


23.84 


2200 


1.5* 


2* 


u 


J101636+563926 


0.232 


108 


24.91 


< 150 


240 


890 


CD? 


J 10573 1+405646 


0.008 


47 


21.59 


1250 


0.5* 


0.1* 


u 


Jl 15727+431806 


0.229 


256 


25.25 


< 150 


630 


2300 


CJ? 


J132513+395552 


0.074 


56 


23.69 


1900 


10 


14 


CSO? 


J 134035+4448 17 


0.065 


82 


23.89 


2300 


3.3 


4.1 


CJ? 


J155927+533054 


0.178 


182 


24.67 


< 150 


1500 


4500 


CSO? 


Additional sources 
















J07 1509+452555 




74 




3800 


1.3 




CD? 


J080454+433537 




360 




1500 


24 




CD? 


J134158+541524 




125 




< 150 


530 




CD? 



only. One source (CORALZ J140942+360416) was too weak at 
5.0 GHz for self-calibration, and we were only able to produce 
a map at 1.6 GHz. 

From these observations we have derived the largest angular 
size, 9, of each source, from which the (projected) largest 
linear size (LLS) can be calculated. Six sources (CORALZ 
J073934+495438, CORALZ Jl 15000+552821, CORALZ 
J140051+521606, CORALZ J143521+505122, CORALZ 
J093609+331308, and CORALZ J 10573 1+405646) were 
(almost) unresolved by our observations, and we have estimated 
6 using the deconvolved major axis. Sizes determined this way 
are less certain and are therefore indicated with a star. The 
results are given in Table Q] with in Col. 1 the source IAU 
name, in Cols. 2-5 the redshift, the flux density at 1.4 GHz, the 
luminosity density at 5.0 GHz, and the peak frequency (Snellen 
et al. 120041 1. in Col. 6 the largest angular size 6, in Col. 7 the 
projected largest linear size, in Col. 8 a possible morphological 
classification (see Sect. 13.21 for a description of the different 
classes). 



4.2. Radio morphology 

For 13 out of the 18 sources from the CORALZ core sample, 
a morphological classification was possible. The majority of 
the classified CORALZ sources show a CSO or CD morphol- 
ogy (10/13, 77%). Only one source has been tentatively clas- 
sified as a core-jet source. These results are in agreement with 
Stanghellini et al. ( fT997l > and Snellen et al. d2000al) . who find 
that GPS sources classified as core-jet are usually identified with 



quasars, whereas GPS sources identified with galaxies tend to 
exhibit CSO morphologies. 

4.3. Spectral turnover and angular size 

Following Snellen et al. (2000b) we searched for correlations be- 
tween the angular sizes and spectral turnovers of the sources in 
the CORALZ core sample, and compared the results with those 
presented in the literature. Early measurements of the angular 
sizes of GPS and CSS sources using VLBI strongly suggested 
that their spectral turnovers are caused by synchrotron self- 
absorption JsSA, Jones, O'Dell & Stein 119741 H odges, Mutel 
& Phillips IT9841 Mutel, Hodg es & P hillips [19851 . This was fur- 
ther advocated by Fanti et al. ( 1990), who showed that the peak 
frequency and the largest angular size of CSS radio sources are 
strongly anticorrelated, as expected for SSA. For a homogeneous 
synchrotron self-absorbed radio source, the peak frequency, v p , 
is given approximately by (Kellermann & Paulinv-Tofh |198lV 

v p ~ 8 B 1/5 Sl /5 <T 4/5 (1 + z) 1/5 GHz, (2) 

where B is the magnetic field in gauss, S v the flux density at the 
spectral peak in Jy, 8 the largest angular size in mas, and z the 
redshift of the source. In contrast, Bicknell et al. ( 1 19971 1 argue 
that models where the radio source is surrounded by ionized gas 
can also reproduce the observed relation, and the turnover could 
also be caused by free-free absorption (FFA). 

Now that we have measured sizes and peak frequencies of 
the CORALZ core sample, at significantly lower redshifts and 
radio luminosities than the literature samples, it is interesting 
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Fig. 1. (left) Largest angular sizes and peak frequencies of the CORALZ core sample (filled cicles) and the bright PHJ (Snellen et 
al. 2002), Stanghellini (1998), and Fanti (1990) samples (open circles). Circles with arrows denote upper limits. The line is a linear 
least-squares fit (in log-log space) to the bright literature samples, (right) To correct for the flux density-dependence predicted by 
SSA, we divide the peak frequency by S^ 5 , and show that the linear sizes and peak frequencies of the CORALZ sources are now 
consistent with the linear least-square fit to the bright samples. 



to see whether weak radio sources obey the same anticorrela- 
tion between 6 and v p . We have collected angular sizes and peak 
frequencies of several flux density-selected samples from the lit- 
erature for comparison. These included the Parkes half-Jansky 
(PHJ) sample of GPS gala xies (S nellen et al. 120021 with an- 
gular sizes from Liu et al. [2007), the GPS galaxies from the 
Stanghellini ( 1998) sample, and the CSS galaxies from the Fanti 
d 1 990 ) sample. The angular sizes for the sources in the latter two 
samples come from O'Dea ( I1998I ). 

Fig- ffl (left) shows that the peak frequency and linear size of 
the CORALZ core sample (filled circles) are also anticorrelated, 
but they fall significantly below the relation for the brighter and 
more distant literature samples. However, this is actually ex- 
pected from SSA theory. As can be seen in equation @, the peak 
frequency of a synchrotron self-absorbed radio source scales 
with the flux density as: v p oc S 2 / 5 , and therefore the low flux 
density CORALZ sample should have smaller angular sizes for 
similar peak frequencies. In Fig.\l\(right) we have corrected for 
this effect, and show that the linear sizes and peak frequencies of 
the CORALZ core sample are now consistent with the linear fit 
to the bright literature samples, perfectly in line with synchrotron 
self-absorption theory. 

In Fig.\l\(right) we assume that GPS and CSS radio sources 
obey equation (O, although this formula in principle only ap- 
plies to homogeneous synchrotron self-absorbed sources. From 
VLBI observations it is clear that GPS/CSS sources are not ho- 
mogeneous, so the fact that the correlation between peak fre- 
quency and overall angular size is so tight, suggests that the ratio 
of component to overall size does not change significantly over 
the lifetime of GPS/CSS radio sources. This fixed ratio of com- 
ponent to overall size is in agreement with VLBI observations 
and can be explained as self-similar evolution (e.g. Snellen et al. 
2000b). To determine how the data constrain the correlation be- 
tween turnover frequency, peak flux density and angular size to 
be in agreement with SSA and self-similar evolution, we fit the 
data using the fitting function: 

e = cS a v ^(l+z) ll \ (3) 

where a, b, and c are parameters that can be varied in order to 
best fit the data. For a = 1/2, b = -5/4, and c = 8 5/4 B 1/4 , this 
expression reduces to equation (f2j). In Fig. [2] we show a con- 
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Fig. 2. Contour map of the goodness-of-fit parameter A^ 2 for the 
combined sample of young radio sources (the CORALZ core 
sample, PHJ (Snellen et al. 2002), Stanghellini (1998), Fanti 
(1990)). The contours are drawn at the 68.3%, 95.4% and 99.7% 
levels of a A^-distribution with two degrees of freedom. The 
location of the parameters expected for SSA (equation (f2]i) is 
indicated by the cross. 



tour map of the goodness-of-fit parameter A% 2 , resulting from 
fitting the data to equation (01 with a range of fixed values for 
a and b. We again used the CORALZ core sample, and the PHJ 
(Snellen et al. 2002), Stanghellini (1998), and Fanti (1990) sam- 
ples. Note that the five larger sources from the CORALZ core 
sample, which have upper limits to the turnover frequency (see 
Fig.[TJ, are not included in this analysis. When turnover frequen- 
cies and peak flux densities for these become available, we ex- 
pect to significantly improve the constraints on the parameters in 
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equation ([3). The contours are drawn at the 68.3%, 95.4% and 
99.7% levels of a A^ 2 -distribution with two degrees of freedom. 
The values of the parameters a and b expected for a homoge- 
neous synchrotron self-absorbed source (equation d2)) are indi- 
cated by the cross. From Fig.[2]we conclude that the data are in 
remarkably good agreement with the homogeneous SSA source 
model, with a = 0.45 + 0.25 and b = -0.9 + 0.3. This confirms 
that the spectral turnovers in GPS and CSS spectra are generally 
caused by SSA, and also favours self-similar evolution models 
of young radio sources. 
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Fig. 3. Rest-frame turnover frequency against linear size for the 
CORALZ core sample (filled cicles) and the PHJ (Snellen et al. 
2002), Stanghellini (1998), and Fanti (1990) samples (open cir- 
cles). The lines are models from Bicknell et al. (1997) for jet en- 
ergy fluxes of 10 46 , 10 45 5 , 10 45 , and 10 427 erg s _1 (dashed, dot- 
ted, solid, and dot-dashed line, respectively), density power-law 
slope of/3-2, and number density at 1 kpc of 10 cirT 3 . 

Can FFA also explain the relations found? Bicknell et al. 
( 1 19971 1 present radio source evolution models, where the AGN 
ionizes the gas surrounding the radio source by shocks and pho- 
toionization. These models predict an anticorrelation between 
the rest-frame peak frequency and the largest linear size of GPS 
and CSS radio sources, consistent with the data available at 
that time (the samples of Fanti et al. (1990) and Stanghellini 
et al. (|1998l l). To see whether the sizes and peak frequencies 
of the CORALZ core sample are also in agreement with this 
model, we have reproduced Fig. 7 (top right) from Bicknell 
et al. (|T997j l in Fig. [3] and added the data from the CORALZ 
core sample. Sources from the CORALZ core sample are repre- 
sented by filled cicles and sources from the PHJ (Snellen et al. 
2002), Stanghellini (1998), and Fanti (1990) samples by open 
circles. The lines are predictions of the radio source models from 
Bicknell et al. (1997) for jet energy fluxes of 10 45 , 10 45 5 , and 
10 46 erg s (solid, dotted, and dashed line, respectively), den- 
sity power-law slope of /3 = 2, and number density at 1 kpc of 
10 cirT 3 . The fact that the lines fit the literature samples fairly 
well is no surprise, since the model is designed to do so. 

The median radio luminosity of the CORALZ core sample 
is a factor of 630 lower than that of the bright literature sam- 
ples, meaning that the jet energy fluxes should also be a fac- 
tor of 630 lower, assuming a constant fraction of the jet energy 
radiated as synchrotron emission (as expected from analytical 
models, e.g. see appendix Sect. [A}, resulting in a jet energy flux 
of ~ 10 42 7 erg s _1 . This is indicated by the dot-dashed line in 
Fig. [3] All-but-one of the CORALZ core members with mea- 



sured turnover frequency fall below this line, indicating that the 
model is not consistent with these data, or that, for some reason, 
the ratio of jet energy flux to radio luminosity is significantly 
lower for the CORALZ core sample. Of course, it will be pos- 
sible to add extra parameters to the Bicknell model, such that it 
fits the peak frequency and angular size data of GPS and CSS 
sources. However, since the SSA self-similar evolution scenario 
is so straightforward and fits the data well, we believe it is sci- 
entifically more elegant. 

4.4. Expansion velocity and radio power 
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Fig. 4. Expansion velocities of Compact Symmetric Objects in 
units of h~ x c, and radio luminosities at 5 GHz in erg s _1 . The 
upper limit indicated with the cross is our 3<x upper limit to the 
expansion velocity of CORALZ J083 139+460800, the other ar- 
rows indicate Icr upper limits from previous studies. The solid 
line shows the expected relation between expansion velocity and 
radio luminosity, derived from an analytic radio source evolu- 
tion model described in the appendix in Sect. |Aj with arbitrary 
scaling. The dashed line represents this relation scaled down by 
a factor of Vl - 0.9 2 0.44, such that 90% of a sample of ran- 
domly oriented sources that intrinsically would follow the solid 
line, should be located above the dashed line. 

Three sources have been observed twice (March 2000 and 
May 2004) with the same VLBI array at the same frequency, 
allowing investigation of possible structural changes, but only 
one source, CORALZ J083 139+460800, is suitable for deter- 
mination of the expansion velocity (see above). However, we 
found no change in the position (x or y) to within 18 fias, which 
corresponds to a 3<x upper limit to the projected expansion ve- 
locity of v exp < 0.07 h~ l c - 0.10 c. We noted that this upper 
limit is lower than the expansion velocities previously found for 
bright archetypal GPS/CSOs, and we therefore investigated the 
dependence of expansion velocity on radio luminosity. We have 
searched the literature for detections of and upper limits to ex- 
pansion velocities in CSOs. Since the first detections of lobe 
proper motions in 0710+439 (Owsianik & Conway [19981 1 and 
0108+388 (Owsianik, Conway & Polatidis 1998) expansion ve- 
locities have been measured for only a few other sources. In 
Table|2]we summarize these results, and also include upper lim- 
its to expansion velocities for which no significant expansion 
has been measured so far. Col. 1 gives the source IAU name, 
Col. 2 the redshift, Col. 3 the lobe-to-lobe, projected expansion 
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Table 2. Expansion velocities of Compact Symmetric Objects 
presented in the literature. 



IAU Name : v exp No. of epochs Ref. 

li^c 

1 
1 
1 
2 
3 
1 
1 
4 
1 
5 
5 
6 
5 
7 



References: (1) Polatidis & Conway (2003 ); (2) Tay lor et al. J2000t : (3) 
Stangh ellini et al. d2002t; (4) T schager et al. ( 12000b ; (5) Gu gliucci et al. 
(2005 1; (6) Tingay et al. ( 120021 ; (7) Tzioumis et al. ( fT998l >. 



velocity in units of he, converted to the adopted cosmology 
(WMAP5, see Sect. [TJ), Col. 4 the number of observing epochs 
and the years of the first and last observations, and Col. 5 gives 
the reference for the expansion velocity. Note that the upper lim- 
its as presented in Gugliucci et al. (2005 ) had not been corrected 
for cosmological time dilation. The values we present here have 
been corrected and represent rest-frame quantities. 

In Fig. |4] we show the expansion velocities measured in 
GPS/CSOs as a function of radio luminosity. Both the radio 
luminosity and the lobe-to-lobe, projected expansion velocity 
have been corrected for cosmological effects, and represent 
rest-frame quantities. The upper limit indicated with the cross 
is our 3<t upper limit to the expansion velocity of CORALZ 
J083 139+460800, the other arrows indicate lcr upper limits 
from the literature (see Table|2]l. 

Although there is a large scatter, a positive correlation with 
radio luminosity can be seen. A Spearman rank correlation anal- 
ysis shows that, ignoring all measured upper limits to the expan- 
sion velocity, the correlation is significant at a 98% confidence 
level (n = 9, p = 0.717, one-tailed). With the upper limits in- 
cluded in the analysis, the correlation reaches confidence levels 
> 98%, regardless of the precise way we treat these upper limits. 
Furthermore it seems that the scatter can mostly be explained in 
terms of projection effects, since the majority of the sources lie 
between the solid and dashed lines, comprising a range of 90% 
probability. 

Note that we do expect a correlation between expansion ve- 
locity and radio luminosity from radio source evolution the- 
ory. In the appendix in Sect. [A] we present an analytic model 
of a young, expanding radio source with a constant jet power 
in a power-law density environment. This model predicts that 

1 /3 

Vexp ^adio- The solid line represents such a relation, with ar- 
bitrary scaling. The dashed line represents this relation scaled 
down by a factor of Vl - 0.9 2 - 0.44, such that 90% of a sam- 
ple of randomly oriented sources that intrinsically would follow 
the solid line, should be located above the dashed line. 



5. Conclusions 

In this paper we present VLBI and MERLIN observations of the 
CORALZ sample, a nearby sample of 25 compact (6 < 2") radio 
sources, of which 18 form a 95% statistically complete sample in 
the redshift range 0.005 < z < 0.16. We have measured the sizes 
of the sources, and performed a morphological classification. 

The CORALZ core sample follows the well established rela- 
tion between radio spectral peak frequency and angular size, but 
with significantly smaller sizes at any particular peak frequency, 
compared to more powerful and more distant GPS/CSS sources. 
We show that this is exactly as expected from synchrotron self- 
absorption theory, in which the angular size is also proportional 
to the square root of the peak flux density. By least-squares fit- 
ting of the combined sample of CORALZ and bright literature 
samples, we recover the dependencies of angular size on peak 
frequency and peak flux density as expected for synchrotron self- 
absorbed radio sources evolving in a self-similar way. Current 
models that invoke FFA to explain the spectral turnovers in GPS 
and CSS sources (Bicknell et al. 1997) can not explain the rela- 
tively small angular sizes found for the CORALZ core sample. 
We therefore conclude that, although FFA may play a role in 
some sources, the turnovers must be caused by SSA for the ma- 
jority of GPS and CSS radio sources. 

In addition, based on two epochs of global VLBI, we 
have derived a strong upper limit to the expansion velocity of 
CORALZ J083 139+460800 (v exp < 0.07 hr x c). We show that, 
in comparison with more powerful young radio sources, the ex- 
pansion velocity is low, indicating that the expansion speed of 
GPS radio sources is correlated to their radio luminosity. We 
show that this is expected from analytic radio source evolution 
modelling, which predicts v exp oc L 1 ^. (appendix Sec. [A}. The 
scatter in the correlation is consistent with being due to projec- 
tion effects. 
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Appendix A: Radio source evolution modelling 

In this section we present analytic models predicting how the 
expansion velocity of young radio sources depends on radio lu- 
minosity L ra dio, source age f, and environment. Our consideration 
is based on the radio source evolution models of Kaiser & Best 
(2007), hereafter KB. If the sources are physically small (less 
than a few kpc in size) and their radio luminosity is measured 
at a frequency where self-absorption is not important, they are 
in the 'synchrotron loss dominated' regime as discussed in KB 
(Sect. 2.3). In this regime we expect the radio luminosity L ra dio 
to be proportional to Q, the jet power, which is assumed to be 
constant with time, and independent of the density of the source 
environment. If we assume that the density p in the source envi- 
ronment follows a power-law distribution with respect to r, the 
distance to the centre, then 

p= P0 L^y, (A.i) 

where po is the density at an arbitrary reference point ao, and /3 
is the slope of the density profile. Note that po and «o are not in- 
dependent, and the model can only depend on their combination 
poo^. Using equation A2 from the appendix of KB, the linear 
size of the source, D, is given by 



D cc 



Q 



1/(5-/3) 



,3/(5-0 



(A.2) 



This means that the expansion speed of the lobe is given by 
v exp = D oc Q m> t (fi-W ( po ^)- 1/(5 -« . (A.3) 
Assuming Q is proportional to L ra dio, we obtain 

^^r^'i^-'^r (a.4) 

If we take the often adopted value f3 — 2 (e.g. Bicknell et al. 
1 19971 1, the time dependence drops out and equation JA.4-b re- 
duces to 

V ex P <* *fcL (P()flo) V3 • (A.5) 

We have used this equation for comparison with the observed 
expansion speeds in Sect. 14.41 Note that in the central part of 
galaxies one might expect the density profile to become less 
steep, even possibly resulting in a flat density profile with fi = 0. 
In this extreme case you would expect v exp oc l}'tt _2 ^ 5 p " 1 ^ 5 , 
so the slope of the correlation between expansion velocity and 
radio power is not extremely sensitive to the value of /3, but the 
additional time dependence may increase the scatter around it, in 
addition to the scatter caused by projection effects and by vari- 
ations in the density profile of the surrounding medium from 
source to source. 
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Appendix B: Observations and radio maps 

Table B.l. Observations of the sources in the CORALZ sample. The first 18 objects form the CORALZ core sample at 0.005 < z < 
0.160, the next seven objects are other nearby radio sources in the CORALZ sample, and the final three objects are other compact 
radio sources located towards probably random foreground galaxies. 
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ED 
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E9 
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1.658 
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-31 
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14 


0.08 
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13 
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ED 
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1.408 
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0.09 
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-47 
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ED 


2000/03/02 
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-20 
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32 
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15 
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EH 


2000/02/14 


1.655 


32 


EVN 


182 


0.39 


27 x 14 


39 


J134158+541524 


EH 


2002/02/19 


1.650 


16 


EVN 


10 


0.13 


25 x 22 


70 



J0733+56 EVN Only 1.6 GHz (2000) 
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J1 037+43 EVNOnly 1.6 GHz (2001) 
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J1718+54 EVNOnly 1.6 GHz (2000) 
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Fig.B.2. Contour maps of the three most compact sources from the CORALZ core sample taken with global VLBI at 5.0 GHz in 
2000 (first column), and 2004 (second column). 
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Fig. B.4. continued. 




Fig.B.5. Contour maps of the additional sources, observed with global VLBI (5.0 GHz), or with the EVN (1.6 GHz). These sources 
have been excluded from the CORALZ sample because the offset between the radio and optical position is large, indicating that 
these radio sources are located towards probably random foreground galaxies. 
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Table B.2. The fitted parameters of the observed components. 



IAU name 


Fi 


g Epoch 


Vobs 


Ax 


Ay 


Flux density Size F 








(GHz) 


(mas) 


(mas) 


(mJy) 


(mas) 


J073328+560541 


IB. 11 


2000/02/14 


1.655 


0.0 


0.0 


248.0±0.9 


29.4x8.4 




ib~ti 


2004/06/05 


1.665 


0.0 


0.0 


150.0+0.7 


7.1x3.8 










-6.7 


15.3 


60.1+0.5 


5.5x2.2 










-4.3 


28.9 


28.6+0.7 


5.9x3.7 




eh 


2004/05/24 


4.993 


0.0 


0.0 


58.9+0.1 


5.1x3.4 










-6.4 


15.5 


28.1+0.1 


5.3x1.8 










-5.8 


26.8 


14.8+0.1 


6.3x4.8 










7.5 


-18.6 


3.1+0.1 


3.1x0.0 


J073934+495438 


15721 


2000/03/02 


4.987 


0.0 


0.0 


30.0+0.1 


1.8x0.4 




15721 


2004/05/24 


4.993 


0.0 


0.0 


36.7+0.1 


2.0x1.0 


J083 139+460800 


E2 


2000/03/02 


4.987 


0.0 


0.0 


54.1+0.2 


1.0x0.8 










-0.9 


-4.3 


19.9+0.2 


0.8x0.1 




15721 


2004/05/24 


4.993 


0.0 


0.0 


65.0+0.2 


1.1x0.7 










-0.9 


-4.3 


35.5+0.2 


0.9x0.6 










0.8 


4.2 


2.3+0.2 


0.6x0.0 


J083637+440109 


15731 


2000/11/13 


1.658 


0.0 


0.0 


60.4+0.8 


385.9x295.2 










546.3 


-5.4 


46.7+0.3 


212.6x64.4 


J090615+463618 


15711 


2001/05/27 


1.659 


0.0 


0.0 


160.7+0.5 


9.5x3.3 










-1.3 


-22.5 


65.1+0.5 


11.1x4.2 




15711 


2004/06/05 


1.665 


0.0 


0.0 


160.7+0.5 


9.5x3.3 










-1.3 


-22.5 


65.1+0.5 


11.1x4.2 




ib~ti 


2004/05/24 


4.993 


0.0 


0.0 


105.0+0.3 


1.1x0.3 










0.4 


6.4 


8.5+0.5 


2.5x1.4 










2.6 


-5.5 


10.3+0.6 


2.7x1.6 










-0.4 


-7.4 


2.2+0.4 


2.8x0.2 










3.4 


-24.6 


6.4+0.5 


1.7x1.3 


J102618+454229 


I57T1 


2001/05/28 


1.659 


0.0 


0.0 


111.0+0.3 


17.9x2.9 




ib~ti 


2004/06/05 


1.665 


0.0 


0.0 


59.0+0.2 


3.5x1.8 










5.7 


-3.3 


9.6+0.2 


1.5x0.0 










-6.4 


4.1 


13.8+0.3 


5.8x3.9 










10.4 


-5.0 


22.8+0.4 


8.7x3.4 




IBTD 


2004/05/24 


4.993 


0.0 


0.0 


33.7+0.1 


0.8x0.3 










5.3 


-1.7 


21.7+0.3 


3.4x1.8 










-4.8 


2.2 


26.0+0.3 


3.0x1.8 










8.8 


-3.3 


5.5+0.2 


1.4x1.1 










1.7 


-0.6 


2.2+0.1 


0.7x0.5 










-7.1 


3.6 


8.9+0.2 


2.6x1.0 


J103719+433515 


ib~ti 


2001/05/27 


1.659 


0.0 


0.0 


72.0+0.3 


13.4x8.8 




ib~ti 


2004/06/05 


1.665 


0.0 


0.0 


24.7+0.2 


3.0x1.6 










-3.0 


2.9 


11.7+0.2 


4.5x0.5 










4.1 


0.7 


24.3+0.3 


6.2x2.1 










-7.2 


6.0 


25.1+0.3 


6.2x4.8 










9.3 


-2.4 


10.4+0.3 


5.6x3.1 




ib~ti 


2004/05/24 


4.993 


0.0 


0.0 


21.0+0.1 


4.2x1.8 










5.5 


-3.2 


13.6+0.1 


1.2x0.0 










-10.3 


6.3 


12.9+0.1 


7.1x5.2 










11.4 


-8.8 


5.1+0.1 


8.6x0.0 


Jl 15000+552821 


15731 


2001/05/28 


1.659 


0.0 


0.0 


111.9+0.3 


7.0x4.1 










40.6 


0.0 


5.8+0.7 


33.5x28.9 


J120902+411559 


IbTTI 


2001/05/27 


1.659 


0.0 


0.0 


128.5+0.2 


8.1x0.0 




I57T1 


2004/06/05 


1.665 


0.0 


0.0 


56.3+0.1 


2.6x0.0 










-2.7 


0.5 


58.9+0.2 


12.2x2.2 










-8.6 


1.0 


25.1+0.2 


4.8x0.0 










4.5 


-2.0 


33.3+0.2 


5.3x3.5 




I57T1 


2004/05/24 


4.993 


0.0 


0.0 


19.7+0.1 


0.6x0.0 










-3.0 


0.2 


7.6+0.2 


1.5x0.8 










-4.5 


0.8 


9.3+0.2 


1.7x1.2 










1.4 


-0.2 


39.0+0.3 


4.2x1.0 


J131739+411545 


1531 


2000/03/02 


4.987 


0.0 


0.0 


126.8+0.6 


1.8x1.2 










1.6 


2.1 


44.0+0.7 


2.8x0.9 



PA 



164 
156 
153 
107 
147 
147 
80 
83 
142 
140 
52 
106 
45 
78 
163 
76 
1 



1 
1 
1 
1 
1 
1 

1 
1 

1 
1 
1 
^ 
1 
1 
1 
1 
1 
1 
1 
1 
1 

c 
1 

^ 
1 

1 
1 
1 
1 

1 

10: 
82 
19 
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Table B.2. The fitted parameters of the observed components, continued. 
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IAU name 


Fi 


g Epoch 


Vobs 


Ax 


Ay 


Flux density Size PA 










(GHz 


) (mas) 


(mas) (mJy) 


(mas) 






IB.2I 


2004/05/24 


4.993 


0.0 


0.0 


183.8+1.3 


2.2x1.9 


49 










2.7 


2.8 


23.3+1.1 


2.3x1.4 


155 


J14005 1+521606 


1531 


2000/12/26 


1.408 


0.0 


0.0 


170.0+0.5 


182.3x95.9 


102 


J140942+360416 


EH 


2001/05/28 


1.659 


0.0 


0.0 


104.1+0.2 


10.7x4.9 


159 










-0.5 


22.1 


22.4+0.3 


22.6x5.4 


16 




ib~ti 


2004/06/05 


1.665 


0.0 


0.0 


111.1+0.3 


13.0x5.8 


169 










1.4 


26.9 


10.3+0.2 


6.1x4.4 


148 


J143521+505122 


1531 


2000/12/26 


1.408 


0.0 


0.0 


75.2+0.6 


161.0x9.3 


69 










128.8 


-46.0 


64.7+0.6 


181.5x46.6 


126 


J150805+342323 


IB3I 


2002/02/19 


1.650 


0.0 


0.0 


4.0+0.2 


7.8x0.0 


72 










-97.5 


-144.4 


6.1+0.4 


42.2x30.5 


47 


J160246+524358 


ib~ti 


2001/05/27 


1.659 


0.0 


0.0 


26.0+0.1 


0.0x0.0 













13.8 


-2.9 


200.3+0.7 


70.4x22.7 


100 










-73.0 


16.2 


222.0+0.8 


64.0x32.8 


106 




ib~ti 


2004/06/05 


1.665 


0.0 


0.0 


41.9+0.4 


3.7x0.0 


145 










19.1 


-7.8 


20.5+0.6 


6.1x2.8 


49 










-23.5 


6.2 


7.4+0.6 


5.5x3.6 


71 










-36.2 


-3.3 


22.1+2.4 


27.3x13.4 


85 










-63.6 


8.8 


26.8+1.5 


15.9x8.3 


132 










-104.6 


43.2 


2.2+0.6 


6.3x3.5 


11 










65.7 


-0.2 


2.8+0.5 


4.9x3.6 


30 




ED 


2004/05/24 


4.993 


0.0 


0.0 


37.7+0.1 


4.1x3.0 


75 










19.1 


-8.7 


23.1+0.2 


8.5x2.6 


105 










-17.5 


8.9 


5.1+0.2 


7.6x2.9 


74 










-37.5 


5.5 


4.8+0.2 


13.0x2.8 


92 


J161 148+404020 


1531 


2000/11/25 


1.658 


0.0 


0.0 


227.2+0.4 


253.2x136.7 


43 










292.8 


-15.9 


214.3+0.5 


304.5x222.0 


22 










-591.8 


-88.8 


22.0+0.6 


329.3x263.1 


134 


J171854+544148 


ED 


2000/02/14 


1.655 


0.0 


0.0 


98.9+0.8 


5.3x2.7 













10.9 


-50.0 


151.4+1.8 


19.2x5.5 


169 




i5~n 


2004/06/05 


1.665 


0.0 


0.0 


193.6+0.5 


6.7x3.0 


9 










9.9 


-45.2 


167.9+0.7 


17.3x2.9 


180 




i5~n 


2004/05/24 


4.993 


0.0 


0.0 


61.5+0.2 


4.7x2.1 


150 










10.2 


-47.4 


26.0+0.2 


5.0x3.2 


115 










9.9 


-36.8 


13.6+0.2 


3.2x1.5 


155 










0.9 


-19.8 


3.8+0.3 


9.1x0.0 


98 










0.4 


20.0 


2.3+0.2 


5.1x0.0 


90 


J093609+331308 


1531 


2000/03/02 


4.987 


0.0 


0.0 


51.7+0.3 


1.5x0.3 


142 


J101636+563926 


1531 


2000/11/09 


1.658 


0.0 


0.0 


57.2+0.2 


64.4x26.5 


154 










136.9 


-196.6 


39.2+0.1 


66.5x10.7 


137 


J 10573 1+405646 


1531 


2000/03/02 


4.987 


0.0 


0.0 


14.8+0.1 


0.5x0.2 


1 


J115727+431806 


1531 


2001/01/02 


1.408 


0.0 


0.0 


43.2+0.2 


55.5x44.2 


76 










-185.3 


604.5 


59.4+0.2 


215.0x131.0 


165 










134.8 


188.5 


30.4+0.2 


100.8x55.2 


10 










41.4 


423.5 


58.4+0.3 


232.8x187.6 


155 


J132513+395552 


1531 


2000/03/02 


4.987 


0.0 


0.0 


16.2+0.1 


0.4x0.0 


176 










2.2 


5.4 


3.0+0.2 


1.0x0.8 


86 










-2.0 


-4.2 


5.6+0.7 


6.1x2.0 


15 


J 134035 +4448 17 


1531 


2000/03/02 


4.987 


0.0 


0.0 


78.3+0.2 


0.2x0.2 


13 










0.7 


3.2 


5.7+0.3 


1.8x0.7 


9 


J155927+533054 


1531 


2000/12/18 


1.408 


0.0 


0.0 


83.0+0.8 


936.4x390.6 


33 










560.2 


1342.1 


63.2+0.5 


555.1x367.6 


4 


J07 1509+452555 


IB. 51 


2000/03/02 


4.987 


0.0 


0.0 


39.1+0.2 


0.6x0.2 


119 










0.5 


-1.2 


11.3+0.2 


0.8x0.1 


162 


J080454+433537 


1531 


2000/02/14 


1.655 


0.0 


0.0 


194.1+0.4 


6.7x3.7 


70 










-24.1 


-1.8 


71.3+0.4 


15.5x3.2 


63 


J134158+541524 


1531 


2002/02/19 


1.650 


0.0 


0.0 


66.3+0.8 


67.9x40.9 


50 










-270.9 


467.1 


10.9+0.4 


33.8x16.3 


108 



